AD-A153  942  VARIA8ILITV  OF  MEASURES  OF  WEAPONS  EFFECTIVENESS  VOLUME 
2  APPLICATION  TO.  .  (U)  FLORIDA  UNIV  GAINESVILLE  DEPT  OF 
INDUSTRIAL  AND  SVSTEMS  ENGIN.  .  B  D  SIVAZLIAN  FEB  85 
UNCLASSIFIED  AFATL-TR-84-92-V0L-2  FB8635-83-C-B282  F/G  12/1 


r 

END  1 

1/ 

NL 


AD-A153  942 


AFATL-TR-84-92 


Variability  of  Measures  of 
Weapons  Effectiveness 

Volume  II:  Application  to  Blast 
Sensitive  Targets  in  the  Absence 
of  Delivery  Error 


BD  Sivazlian 


UNIVERSITY  OF  FLORIDA 

DEPARTMENT  OF  INDUSTRIAL  AND  SYSTEMS  ENGINEERING 
GAINESVILLE.  FLORIDA  32611 


FEBRUARY  1985 

FINAL  REPORT  FOR  PERIOD  MAY  1983  •  JANUARY  1985 


DT1C 

SELECTEW| 
MAY  2  1 1985  ^  9 

B 


Approved  for  public  release;  distribution  unlimited 


O 


P  Air  Force  Armament  Laboratory 

°  AM  FORCE  SYSTEMS  COMMANO*  UNITED  STATES  AIR  FORCE* EGLIN  AIR  FORCE  BASE,  FLORIDA 


V 


Ji 


UNCLASSIFIED _ _ _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


l»  REPOR7  SI  LUflITY  CLASSIFICA1  ION 

UNCLASSIFIED 


2>  SECURITY  CLASSIFICATION  AUTHORITY 

N/A 


2b  DECLASSlFlCATION/OOWNGRAOING  SCHEDULE 

N/A _ 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

N/A 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


6*  NAME  OF  PERFORMING  ORGANIZATION 

University  of  Florida 

6b.  OFFICE  SYMBOL 

I If  applicable ) 

N/A 

6c  AODRESS  Kifv  State  and  /■If  Code l 

Department  of  Industrial  and  Systems 
Engineering 

Gainesville,  Florida  32611 

NAME  OF  FUNOING/SPONSORING 
ORGANIZATION 

Analy  ^  Strat  Def  Div 

8b.  OFFICE  SYMBOL 

Of  applicable ) 

DLY 

Be  ADDRESS  ,Citv.  Stale  and  /.If  Code ) 

Air  Force  Armament  Laboratory 
Eglin  Air  Force  Base,  Florida 

32542 

Iv'Wii 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  Public  Release 
Distribution  Unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUMBERISI 

AFATL-TR-84-92,  Volume  II 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Weapon  Evaluation  Branch  (DLYW) 
Analysis  Division  _ 


7b.  AODRESS  (City.  State  and  /IV  Code! 

Air  Force  Armament  Laboratory 
Eglin  Air  Force  Base,  Florida  32542 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

c 

Contract  No.  F08635-83-0202 


10  SOURCE  OF  FUNDING  NOS 


TASK  j 

WORK  UNIT 

NO 

NO 

25 

03 

12.  PERSONAL  AUThOR(S) 

SiY-^lian 

fl3*  type  of  report 


13b.  TIME  COVERED 

f Ron  May  to  Jan  f 


flZoATEOFREPORTWVrTMo  .  Day >  I  15.  PAGE  COUNT 


16  supplementary  notation 

Availability  of  this  report  is  specified  on  verso  of  front  cover. 


17 

COS  AT  1  COOES 

18.  SUBJECT  TERMS  1 Continue  on  reverse  if  necessary  and  identify  by  block  number i 

field 

GROUP 

SUB.  GR. 

1 1 

_ QJ _ 

Blast  damage  function,  statistical  estimation 

19.  ABSTRACT  iContmue  on  revert*  if  necessary  and  identify  by  block  numberi 

f  The  problem  of  computing  the  uncertainty  associated  with  the  probability 

of  kill  Pj^  due  to  blast  in  the  absence  of  weapon  fragmentation  and  aiming 

error  is  considered.  The  damage  function  due  to  blast  is  assumed  to  be  a  piecewise 
linear  function  of  the  distance  between  the  blast  point  and  the  target  point.  Under 
tiie  assumption  that  the  input  parameters  are  uniformly  distributed,  explicit 
expressions  are  obtained  for  the  expectation  and  variance  of  P]^_ 


20  distribution  availability  of  abstract 
UNCLASSIFIED  unlimited  2  SAME  AS  RPT  C3  OTIC  USERS  U 


122*  NAME  OF  RESPONSIBLE  INDIVIDUAL 

■'-'a  ■  Daniel  Mclnnis _ 

0D  FORM  1473,  83  APR 


21  ABSTRACT  security  CLASSIFICATION 

Unclassified 

22b  TELEPHONE  NUMBER  |22c  OFFICE  SYMBOL 

/ Include  \  rea  Code) 


(904)  882-4455 

EDITION  OF  1  JAN  73  IS  OBSOLETE 


AFATL./DLYW 


SECURITY  CLASSIF  ICA  T  ON  OF  This  PAGE 


UNCLASSIFIED _ 

SeCUHITY  CLASSIFICATION  OF  THIS  PAGE _ 

11.  TITLE  (Concluded) 

Application  to  Blast  Sensitive  Targets  in  the  Absence  of  Delivery  Error 


PREFACE 


Thfs  report  describes  work  done  In  the  summer  of  1983  by  Or  B.  D. 
Slvazllan,  Department  of  Industrial  and  Systems  Engineering,  the  University  of 
Florida,  Gainesville,  Florida  32611  under  Contract  No.  F08635-83-C-0202  with 
the  Air  Force  Armament  Laboratory  (AFATL),  Armament  Division,  Eglin  Air  Force 
Base,  Florida  32542.  The  program  manager  was  Mr  Daniel  A.  Mclnnis  (DLYW). 

The  work  was  initiated  under  a  1982  USAF-SCEEE  Summer  Faculty  Research 
Program  sponsored  by  the  Air  Force  Office  of  Scientific  Research  conducted  by 
the  Southeastern  Center  for  Electrical  Engineering  Education  (SCEEE)  under 
Contract  No.  F49620-82-C-0035. 

This  work  addresses  Itself  to  the  problem  of  computing  the  uncertainty 
associated  with  the  probability  of  kill  Pkb  due  to  blast  alone  in  the  absence 
of  weapon  fragmentation  and  aiming  error.  The  assumption  is  made  that  (1)  Pkb 
is  unity  between  the  center  of  the  blast  and  a  distance  A  from  the  center,  (2) 
pkb  ’s  ne9H9*b1®  beyond  a  certain  distance  B  from  the  center  (B  >  A),  (3)  P^ 
decreases  linearly  between  A  and  B,  and  (4)  A  and/or  B  are  random  variables 
uniformly  distributed  over  given  ranges.  Explicit  expressions  for  E[Pkb]  and 
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E[P£bl  are  derived  for  various  weapon  target  locations.  (E[*J,  as  usual,  is 
the  expectation  operator.) 

The  author  has  benefited  from  helpful  discussions  with  several  people. 
Particular  thanks  are  due  to  Mr  Jerry  Bass,  Mr  Daniel  Mclnnis  and  Mr  Charles 
Reynolds,  all  from  DLYW,  who  have  read  the  report  and  have  contributed  to  it 
through  helpful  comments. 

The  report  is  the  second  of  a  series  dealing  with  the  uncertainty 
associated  with  various  weapon  effectiveness  indices  and  details  methodologies 
and  techniques  used  in  computing  such  uncertainties  in  the  presence  of  error 
in  the  input  parameters. 
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SECTION  I 


INTRODUCTION 


When  for  an  exploding  weapon  the  effect  of  fragmentation  Is  negligible,  a 
target  may  be  damaged  or  destroyed  due  to  the  prevailing  blast.  If  the  target 
is  close  to  the  location  of  the  weapon,  the  probability  that  the  target  will 
be  destroyed  is  unity.  The  effect  of  the  blast  prevails  highly  up  to  a 
certain  distance  from  the  center  of  the  exploding  weapon,  and  after  some 
point,  it  decreases.  The  damage  becomes  negligible  or  non-existent  beyond  a 
certain  distance.  It  is  customary  to  describe  the  probability  of  kill  P^  due 
to  blast  as  a  function  of  the  distance  between  the  center  of  the  weapon  and 
the  center  of  the  target  by  a  graph  as  shown  in  Figure  1.  Analytically,  Pkb 
as  a  function  of  distance  R  has  the  form 

1  0  <  R  <  A 

1  -  A  <  R  <  8  (1) 

B-A 

0  B  <  R  <  ~ 


Pkb 

1 

0 

FIGURE  1.  Probability  of  Kill  due  to  Blast  as  a  Function 
of  the  Distance  Between  the  Center  of  the  Weapon 
and  the  Center  of  the  Target. 
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A  Is  the  maximum  distance  at  which  irreparable  damage  is  experienced,  and  B  is 
the  minimum  distance  at  which  no  damage  occurs.  It  is  assumed  that  Pkb 
decreases  linearly  between  the  distances  A  and  B. 

The  distances  A  and  B  are  obtained  experimentally  for  a  given 
weapon/target  combination  and  are  subject  to  measurement  error.  Thus,  A  and  B 
are,  in  general,  random  variables.  The  nature  of  the  distribution  of  A  and  B 
will  depend  on  several  factors  such  as  the  weapon,  the  target,  the  atmospheric 
conditions,  etc.  In  the  present  work  it  shall  simply  be  assumed  that  if  A  is 

a  random  variable,  it  is  uniformly  distributed  in  the  interval  Aj  <  A  <  A?. 

Similarly,  if  B  is  a  random  variable.  It  is  uniformly  distributed  in  the 
interval  Bj  <  B  <  B^.  The  assumption  is  further  made  that  the  four  quantities 

Aj,  A2,  B^,  and  B^  can  be  determined  and  that  A^  <  A^  <  B^  <  B2. 

Under  these  assumptions  and  conditions,  the  questions  to  be  answered 

consist  in  estimating  the  mean  and  variance  of  Pkb,  the  probability  of  kill 

due  to  blast.  It  is  evident  that  Pkb  is  a  random  variable  and  that  its 

estimated  values  will  depend  on  the  relative  distance  R  between  the  center  of 

the  target  and  the  center  of  the  weapon.  Hence,  as  may  be  expected,  E[Pkb] 

2 

and  E[Pkb]  will  assume  different  mathematical  forms  depending  on  the  distance 
R.  £[•],  as  usual,  represents  the  expectation  operator.  The  variance  of  Pkb, 
Var[Pkb]  can  then  be  obtained  from  the  well  known  relation 

Var£'W  ■  «pkbi  -  W'W!2 

The  estimates  of  Pkb  will  be  investigated  for  the  following  three 
situations: 

1.  A  is  a  random  variable  and  B  is  fixed. 

2.  A  is  fixed  and  B  is  a  random  variable. 

3.  A  and  B  are  random  variables. 


SECTION  II 


COMPUTATION  OF  E[Pkfa]  AND  E[P^b3 

1.  A  Is  a  random  variable  and  8  Is  fixed 

The  probability  density  function  of  A  Is  assumed  to  be 

0  0  <  x  <  Aj 

fA<*>  ■  |  N  < »  <  h  <2> 

0  A2  <  x  <  - 

Further,  the  distance  B  is  fixed  and  known.  Thus  Aj,  A2,  and  B  are  input 
parameters.  The  center  of  the  target  may  be  located  such  that  its  distance  R 
from  the  center  of  the  blast  satisfies  any  of  the  following  inequalities: 

0  <  R  <  Aj,  Aj  <  R  <  A2,  A2  <  R  <  B,  and  B  <  R  <  »  . 

Each  of  these  four  cases  Is  considered  separately,  and  expressions  for  E[Pkb] 

2 

and  E[Pfcb]  are  computed. 

Case  I:  0  <  R  <  Aj 

From  relation  (1)  it  Is  evident  that  for  all  R  -  lying  between  0  and  the 
ml nlmum  value  of  A,  pkb  =  1,  It  follows  that 

E[Pkb]  =  1  and  E[P^b]  =  1  . 
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Case  2:  A^  <  R  <  A^ 


Figure  2.  Pkb  =  1  -  jr* 

if  Aj  <  A  <  R  <  A, 

pkb  -  1 

if  A1  <  R  <  A  <  A. 

From  relation  (1)  and  from  Figure  2,  it  is  evident  that  pkb  takes  two 
different  values  depending  on  the  relative  position  of  R  with  respect  to 
A^  A  and  A2.  Hence, 


1  -  inr  <  *  <  «  < 

1  Aj  <  R  <  A  <  A2  . 


(3) 


Then  , 


EEPkb]  = 


A 

I2  P 

A1 


kb 


fA  (x)  dx 


R 


JT5P  ttr 


1  dx  +  /  2  1 

D 


dx 


r  pkb  •  M  dx 

Ai 


R  d  *  ^  1  .  /•  2  1 2  1  >jv 

•  /,  &-S#  ^dK  +  '*  1  '*ft 


s  jB  '  R)  f  — - — *•  dx  + 

(A2-A jT  (B-x)2 


a2-r 

A2"A1 


a2-r 


,  (B-R)(R-A,) 

-  jpr^  t"  iB-qr  *  (VR)  1 


ase  3:  A2  <  R  <  B 


For  this  particular  case,  It  may 


be  noted  from  Figure  3  that 


The  probability  density  function  of  B  is  assumed  to  be 


0 

*  \ 

0 


0  <  y  <  Bj 
<  y  <  B2 
B9  <  y  <  •  , 


(10) 


F  < 


The  results  obtained  so  far  for  E[Pkb]  and  E[Pkbl  are  summarized  on  the 

following  page.  Note  that  at  the  values  of  R  =  A^t  R  *  ^2  an<*  R  *  B»  both 

2 

ECPkb]  and  E[P^b]  are  continuous  function  of  R. 


(  1 


E[Pkb]  - 


(A2-Aj5  KB“R)  1n  B-R1  +  (A2"RH 


B  -  R  ,  B"A1 
ApA[  OT 


f  1 


ECP&,]  • 


i  (B-R)(R-A.) 

t  -  Ift-K’jT"  +  (VR)  1 


(B-R ) 

(b-a^Tb-Aj) 


{  0 


2.  A  is  fixed  and  B  is  a  random  variable 


0  <  R  <  A 


1 


Aj  <  R  <  A2 


A2  <  R  <  B 


B  <  R  < 


0  <  R  <  A, 


Al.  <  R  <  A2 


A2  <  R  <  B 


B  <  R  <  - 


(8) 


(9) 


;  .\v 
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Further,  the  distance  A  Is  fixed  and  known.  Thus,  A,  B^»and  B2  are  Input 
parameters.  The  center  of  the  target  may  be  located  such  that  Its  distance  R 
from  the  center  of  the  blast  satisfies  any  of  the  following  Inequalities: 

0  <  R  <  A,  A  <  R  <  Bj,  Bj  <  R  <  B2,  and  B2  <  R  < 

Each  of  these  four  cases  Is  considered  separately,  and  expressions  for  E[Pkbl 
and  E[p£b]  are  computed. 

Case  1:  0  <  R  <  A 

From  relation  (1),  It  Is  evident  that  for  all  R  lying  between  0  and  the 
distance  A,  Pkb  -  1.  Hence, 

E[Pkb]  =  1  and  E[P^b]  *  1  . 

Case  2:  A  <  R  <  B1 


A  R  Bj  B  distance 

Figure  4.  pkb  =  1  *  if  A  <  R  < 

It  may  be  noted  from  the  figure  that  In  this  case 


From  relation  (1)  and  from  Figure  5,  It  Is  evident  that  Pkb  takes  on  two 
different  values  depending  on  the  relative  position  of  R  with  respect  to 
Bj,  B,and  82*  Hence 


Bj  <  B  <  R  <  B2 


Bl  <  R  <  B  <  B2  . 


Then , 


ebw  -  r  pkb  •  fB  (y)  dy 

Bi 


•  ;Bl  0  •  dy  *  V 


Oft 

vv  v*r 


1QT 


Similarly , 


^  -  / 2  p2b  fB  W  * 

B1 


■  /  o2-rVdy  +  /B2  (i-f^)2g-Vdy 

Bj  2  1  R  y  2  1 

*D  R  y  A  (y-A;  B2  *1 

=  17^7  t(B2'R)  '  2(R'A)  ln t A"  ■  (R_A)2  fepr ■ 

.  B--A  (R-A)(B?-R) 

=  rar  KVR)  ‘  2(R'A)  1n  inr +  — • 


Case  4:  B2  <  R  <  « 


Clearly  Pwh  =  0,  hence 


E[Puh3  -  0  and  E[P‘h ]  =  0  . 
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The  results  obtained  so  far  for  ECP^]  and  E[Pk&]  are  summarized  bn  the 

following  page.  Note  that  at  the  values  of  R  =  A,  R  =  Bj»  end  R  =  B2,  both 
2 

ECP^j,]  and  EP^]  remain  continuous  functions  of  the  variable  R. 


EPW  - 


1  R~A  i n  B2‘A 

1  ‘Bpr^1nFp7r 


R-A  ®2-A 

nr 


0  <  R  <  A 
A  <  R  <  B 


1 


B^  R  <  B? 


B2  <  R  < 


(16) 


wk2bi  ■ 


-  2 


R-A  in”2 


B--A 


(r-a  y 


B2-R 

V*i 


B2"B1  10  BrA  +  (B2-A)(B1-A) 


o  R-A  ,  B2”A  (R-a)(B2-R) 

2  B^Bj  ln  R-A  +  (B2-Bx) (B2-A) 


0  <  R  <  A 
A  <  R  <  Bx 

(17) 


Bt  <  R  <  B2 
B2  <  R  <  »  . 


3.  A  and  B  are  random  variables 

It  will  be  assumed  that  the  probability  density  functions  of  A  and  B  are 
given, respectively,  by  (2)  and  (10).  The  Input  parameters  to  the  problem  are 
Aj,  A2,  Bj,  and  B2  where  <  A2  <  <  82.  Again,  here  the  center  of  the 

target  may  be  located  such  that  Its  distance  R  from  the  center  of  the  blast 
satisfies  any  of  the  following  Inequalities: 


0  <  R  <  Ar  A,  <  R  <  »J,  A2  <  R  <  8r  B,  <  R  <  tp  Bj  <  R  < 
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Each  of  these  five  cases  is  considered  separately,  and  expressions  for  E[Pkb] 
2 

and  E[Pkb]  are  computed. 

Case  1:  0  <  R  <  A^ 

Since  for  this  case  Pkb  =  1,  It  follows  that 
E[Pkb]  =  1  and  E[Pkb]  =  1. 

Case  2:  Aj  <  R  <  A2 


It  is  evident  that 


kb 


1  - 


R-A 


nr 


if  Aj  <  A  <  R  <  A2 


If  A1  <  R  <  A  <  A2 


Hence,  from  (4) 


1 


1 


ecp"]*£  w^d‘dy 


®2  1  ^  R-x  1  1 

-  /  Hr  1/  t1  -  rH dx  M 2 1  •  Hr  dxl  rfy 

“2  1  Ax  y  x  "2  ftl  R  ft2  *1 


y-A, 


=  (A2-A1)(B2-B1)  H  ^(y'R)  lnyH  +  (A2_R^  dy 


=  ('va7T(b7-'b7T  U2  (y-R> ln  (y-Ai}  dy 


i 


B, 


-  /’2  (y-R)  In  (y-R)  dy}  ♦  ^ 
B1  VA1 


(18 


Each  of  the  Integrals  on  the  right-hand  side  of  (18)  Is  considered  separately 
For  the  first  Integral  one  obtai ns, after  integrating  by  parts. 


/  2  (y-R)  In  (y-A.)  dy  =  In  (y-A  ) 

B1 


D 

B 


2 

1 


(Bp-R)2  (B,-R)2 

=  - In  (B2-A1)  -  — ^ - In  (B1-A1) 

B  (R-A  )2 

-7  /B  [y  +  Ai  -  2R  +  Tyrjqy  1  d* 

(B?-R)2  (B  -R)2 

•  — ~2 - ln  (Bp-A^  -  — \ - In  (B^Aj) 

i  p  Bp 

-  ^  ^  +  (Ar2R)  y  ♦  (R-A^j  In  (y-Aj)]^ 

(Bp-R)2  (B,-R)2 

=  — — 1n  (W  -  — r — ln  (W 

B  -A 

-  j  (B2-B2)  -  |  (Aj-2R)  (Bp-Bj)  -  |  ( R-Aj ) 2  In  5^-  . 


The  value  of  the  second  Integral  In  (18)  can  be  obtained  from  (19)  by 
substituting  R  to  Aj.  This  yields 


fa  (y-R)  in  (y-R)  dy  = 

Bi 

(Bp-R)2  (B,-R)2 

— ^ In  (Bp-R)  -  — \ 


ln  (B^-R)  - 


B2-B? 


R(Bp-B1) 

- 2 - 


Substituting  (19)  and  (20)  In  (18)  results  In: 


(B  -R)^  (B  -R)2 

EtPkb^  =  fA2-A1) (Bg-Bj)  ^  2  2  ln  *W - ^5 - 1n 


T  ^B2”B1^  -  7  (A1-2R)(B2-B1)  -  {  (R-Aj)2  In 

n2 


B„-A, 


(Bz-R)fc  .  .  (B,-R)^ 


-  [— S~ ln  (b2”r)  ■  — ^2 — ln  (BrR)  - 

R(B--B.) 

+  - +  (A2-R)(B2-B1)  } 


B2  B2 
B2"B1 


which  reduces  to 


(B,-R)2  8,-A,  (B,-R)2  BrAj 


1  UP  n1  \°1  *  i  u  *-n 

ECPki>]  ■  f  5  1 "  sp - 5 —  ' n  BpR 

-7  <R-A1>2  ,n7^t7  (B2-B1)(2S2-ArR)] 


To  compute  E[P^b]  expression  (5)  is  referred  to  and  one  obtains 


(21) 


E[p2bWV  pk2b  dy 


81  A1 


2  1  A?  2 
1  dx  +  J  2  l2 

R 


R  p  Y 

-  J  [J  (1  -  f£)  JT--A- 

«!  A1  A2  A1 

B2  f  (y-R)(R-Aj)  A?-R  i  t 

=  41  ^A2-AiJ(y“AiJ  + 

(R-A,)  B,  „  B  A0-R 


dy 


1'  “2  y-R  ,  n2~ 

/  7TF-  dy  +  T~. 


A  -A 
*2  M1 


dx] 


BrBi 


dy 


'  (A2-A1)(B2-B1)  1  y-Aj  '  7^ 

(R-Aj)  B2  Aj-R  A2-R 

“  (h2'-X1H'b2-b1)  I  1‘  *  pq-1 dy  *  spq- 

•  TXp^TTgpTT  «VBi!  *  (R-Ai>  <y'V 


B?  A?-R 

bJ  +  i^prl 


(R-Aj)^  B2-Aj 

1  ■  Tup'pnr^)  1n^PY* 


(22) 
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Case  3:  A2  <  R  <  Bj 

For  this  case,  Pkb  3  1  ”  *  Thus*  us1ng  (6)  °ne  g*tS 


^  f  (2  dy 


B1  A1 


2  l,v  2  l1  B. 


8?  y'Al 

/  (y-R)  1"  7 tt  dy 


2  ln  2  1  B 


[/  2  (y-R)  In  (y-Aj)  dy  -  /  2  (y-R)  In  (y-A2)  dy  ] 


Using  (19)  one  Immediately  obtains 


E[Pkb]  * 


1 ' k°z_DX J 


j[_2  in  (B2-Aj)  -  -V—  ln  <®r»i> 


B  -A 

{  (BZ-BZ)  -  \  (Aj-BRHBj-Bj)  -  \  (R-Aj)2  '"S7S7 


(B,-R)' 


(B.-R) 


(B,-A  )  -  ^Uj^-ln  (Bj-A2)  -  |  (BZ-BZ) 


2 -  '  '“2  2 


i  (A2-2R)(B2-Bj)  -  \  (R-A2)2  In  rxll 


♦  j  (Bj-BjHAj-Aj)  -  (R-Aj)  ln 


The  expression  for  EtPj^j] 


«*w'*  (2  [i  -  ^  * 


Using  (7)  results  In 


^S*J '  for  f  foffi-v dy  • 


(24) 


The  Integrand  In  (24)  Is  reduced  into  partial  fraction  as  follows 


(v-R 


(y-A2) 


for 


(A. +A?-2R)y-A. A^+Rp 

■'»  (y-vfor- 


=  1  + 


I _ .  t 


where  E  and  F  are  constants  whose  values  are 


E  = 


F  = 


(Ax+A2-2R)  Ap^+R2  (R-A2)2 


t  ~:r 
2  Al 


A  _'A" 
A2  A1 


(A1+A2-2R)  Aj-A^+R2  (R-At)2 


A1'A2 


Al"A2 


So  that 


(v-R) 

{y-T-)(y-\) 


=  l  + 


(R-a2)‘ 

T-_r- 
a2  Al 


y-^ 


(R-Aj)' 

""A"  _A 
2  1 


The  expression  in  (24)  reduces  to 


.  1  8 2  (R-A2)2  1  (R'A1)2  1  , 

E[pkbl  =  ET^7  /.  t1  +  •  yT^  -  'A~Ai  ■  •  yliq-l  dy 


2  1  B 


1 


=  1  + 


(R-A ?r  B  -A 

*  In  6  d 


(R-A,)  B--A. 

1  ln*i-A.  (25) 


(A2-A1TrB7-B7T  (a jpTKspqr ,n*p7  • 
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Case  4:  81  <  R  <  82 

For  this  particular  situation 

<  B  <  R  <  B2 
Bj  <  R  <  B  <  B?  . 

The  expression  for  E[Pkh]  is 


kb 


1  - 


R-A 

8^4 


If 

if 


A„  B„ 

E[pkb]  -  ; 2  t; 2  o 


'1 


'1 


1  .  A  f82  R-x  v 

VB1  {  1  '  1 


VB1 


dy } 


A2'A1 


dx 


Using  expression  (14)  one  immediately  obtains 


A,  B9-R  d  B9-x  , 

E[pkb]  =  j  [  1  n  t r=H  d* 

=  (R"x)  1n  ^2’X>  dx 


1 


-  }  (R-x)  1 n  (R-x)  dx  ]  . 

A, 


Consider  the  first  integral  In  (26)  and  use  integration  by  parts  to 
obtai n 


J  2  (R-x)  In  (B2-x)  dx  =  -  In  (B2-x) 

A1 

(R-A?)2  (R-A.)2 

■ - £■= —  1  n  (B2-A2)  + - gi—  In  (Bg-Aj) 


fA2  (R-x)2  . 


1  a2  (B--R)£ 

+  7  V  [x  +  82  ■ 2R  ■  T^rr]  dx 
A1  z 


(26) 
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(R-Ap)  (R-A  )2 

- 2 -  1n  (82-A2)  +  - T~  1n  (Bg-Aj) 


+  7  [  7~  +  (B2‘2R)  x  +  ^B2-R^2  ln  (Bp-x)  ]|  2 

A1 

(R-A2)2  (R-A.)2 

- - f—  In  (B2-A2)  + - f—\<\  (B2-A1 ) 


+  T  (A2“A?>  +  I  (B2-2R)(A2-A1)  +  \  (B2-R)2  ln  . 


Consider  now  the  second  integral  in  (26).  The  value  of  this  integral  can  be 
obtained  from  (27)  by  substituting  R  for  82.  Thus, 

h  (R-A?)2  (R-A,)2 

/  (R-x)  ln  (R-x)  dx  = - p—  ln  (R-Ap)  +  - yi—  ln  (R-Aj) 


+  {  (A2  -  A2)  -  Y  R  (A2'Al) 


The  expression  for  E[Pkb]  in  (26)  becomes  after  using  the  results  in  (27)  and 
(28) 


E^kb^  “  (A^AjHBg-Bj)  {(82-R ^  ( A2_Ai  ^  -  [  " 

(R-AJ2  ,  o  o 


(R-Az)‘ 


ln  (Bp-Ap) 


+  - r~  ln  (B2_Al)  +  T  (A2  "  Al>  +  7  (82-2R)(A2-A1) 


1  2  ^2~^?  ( R-A^) 

+  2  <VR>  ln  B~Af]  r~1n  (R-A2) 


(R-A.) 


+  - r—  In  (R-Ax)  +  -i(A2-A2)  -^R  (Ap-A^  } 


'^(Bp-Bj)  t  (YA1)(VR)  + 


(R-A 7)c  R?-Ap 

F  ln 


(R'A1  )  1nB2'Al  1  ,n  r,  J  ,  B2"A2  l 
> — 1n  iprp  •  7  <VR) 
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The  expression  for  E[p£b]  is 


E[PkZb3  - 


%  V’siVyt/!!2 11 dyl*^r 


Using  expression  (15)  one  obtains 


.  ,  A,  .  B?-x  (R-x)(B?-R) 

£  T  _  1  f  C  i  f  /  r»  n  \  of  n  \  1  .  ^  ^  


EW  •  r  8p!7  «VR>  -  2(R-‘>  "ET* -Tv^ 


i  ( B^-R ) ( Ap-A. )  o  ^2 

I  '  v;-t—  -  dr  U.  (R-*>  (V*> 


2  1 


'2  “1 


'2  “1  A 


1 


B2“R  A 2  R-x 


2  JZ  (R-x)  In  (R-x)  dx  ]  +  bT-B7  K  &dx}  * 
A1 


2  1  \  2 


Using  (27)  and  (28)  one  obtains 


!  (Bj-RMAj-A,)  2  (R-Aj)' 


Ecp?bi  ■  i-f  b2^- --  -  spr^  I-  — T- 


2  2 


(R-A, )  1971 

+  - In  (Bg-Aj)  +  j  (ApAj)  +  £  (B2-2R)(A2-A1) 

,2 


4(v*>2i"K1  +  ^['^1n(R'A2) 


(r-a2)‘ 


(R-A.) 


+  - J—  1n  (R-Ai)  +  i  (A2'Ai)  -  \  R(VAi)  ] 


b2-r 


B  -A 

*  H»r‘ii  ♦  (r-b2>  ’"ispcjll  • 


Upon  simplification,  the  following  result  is  obtained 
Bo-R 

~  +  TT - X  TTB - B — T  f ( R“A0 ) ^  In  X“T 

2 


b2-a2 


E[pkb3  =  +  (ApA-p'^Y^T  Ur“A2)2  ln  TT 


2  ®?”A1 
(R-A,)Z  In -It1]  • 


Case  5:  Bg  <  R  <  ® 

Here  Pkb  =  0  and  It  follows  that 
E[Pkb]  -  0  and  E[pj|b]  =  0. 

The  results  obtained  so  far  for  E[Pkb]  and  E[Pkb]  when  A  and  B  are  random 
variables  are  summarized  on  the  following  pages.  Note  that  both  E [Pkbl  and 
ECp?b3  are  a  continuous  function  of  R. 


E[Pkb] 


(1)  1 


0  <  R  <  A 


1 


(2) 


1  v  2  '  2 

(A2-A1)'('B2-B1)  t-  2  ln  B^R 


(8,-R)2  Bg-Aj  (B^R)2  B,-A 


1  n 


1  "1 
Bj-R 


Bo-A. 


-  \  (R-Ax)2  ln^pq-  +  \  (B2-B1)(2A2-A1-R)  ] 

{3)  (A2-A1)Tff2'-B1'y  l — 2 — - ? 

B  -A 

♦?<VBi>iy\)  -i<R-Ai>2,nBpf7 


Aj  <  R  <  A2 


(b9-r)2  b,-a,  (B,-R)2  b1-a1 

BrA2 

Bo-A, 


A2  <  R  <  Bx 


(R-A,)2  B,-A,  (R-A,)2  Bg-Aj 


1  r'  2'  i  _  2  2  '""l' 

(4)  ( a2-a"1')Tsj'-b" )"  L— r~1n  R^  r 


1  n 


inf 


B„-A, 


i<VR>2  <  VV<VR> 


1 


Bx  <  R  <  Bg 


(51  0 


B,  <  R  < 


ECPkV 


(1)  1 


Bo-A, 


(R-A.)  u?_r,l 

(2)  1  "  (A2-A1)(B2-B1)  ln 


B„-A, 


( R-A. )  up-n. 

(3)  1  '  (Ag-AjHBg-Bj)  1n  B^ 


(R-A2)' 


(a2-a1)(b2-b1) 


ln 


VA2 

BrA2 


(4) 


b2-r 


(R-Aj)* 


B2‘B1  ’  (A2‘A1){W  1n  ~TZKl 

,2 


BrAl 


(R-A2)‘ 


ln 


B2"A2 


<VA1  (|  -f  )  R-* 


(5)  0 


0  <  R  <  A, 


Aj  <  R  <  A2 


A2<  R  <  Bj 


B1  <  R  < 


B2  <  R  <  - 


(32) 


4.  Example 

Assume  that  an  explosive  charge  detonates  and  that  it  is  required  to 
determine  the  impact  of  the  resulting  blast  on  a  target  located  19  ft  away 
from  the  center  of  the  blast.  The  effect  of  fragmentation  is  neglected.  It 
Is  estimated  that  if  the  target  is  located  at  a  distance  of  15  f t  ^  2  ft,  the 
damage  due  to  blast  is  irreparable.  While  if  the  target  is  located  at  a 
distance  of  21  ft  +_  1  ft,  the  damage  is  negligible.  It  is  required  to 
estimate  the  probability  that  the  target  will  be  damaged. 

Assume  that  A  is  uniformly  distributed  in  the  interval  13  <  A  <  17, 
and  that  B  Is  uniformly  distributed  in  the  Interval  20  <  B  <  22.  Clearly 
Aj-13  ft,  A2«17  ft,  Bj-20  ft,  B2«22  ft,  and  R*19  ft.  Note  that 
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<  R  <  Bj.  This  corresponds  to  the  case  when  A  and  B  are  random 
variables  and  A2  <  R  <  Bj.  Formulas  (31)-3  and  (32)-3  are  pertinent 
and  one  obtains  for  E[P|^]  and  E[Pj^j]  the  following 


ECPvk]  - 


(B,-Rr  B,-A. 


a2-«1md2"d1j 


[-At— in 


B2"^2 


(B ,-Rr  B  -A.  . 

^2—  "’^4(V61>(VA1> 


-?<«V 


22-13 

WAT 


1nf^!T  +  ?  (22-2°)(17-13) 


?  (19-13) ^  in  + 1  (19"l7)  ln  Iferrl 


■F[!,"!-71"T+4-:f,nf  +  7,nt] 

E[Pkb]  =  1  [(4. 5) (.5878)  -  \  (.8473)  +  4  -  18( .2513)  +  2( .5108)  ] 
=  ^  (2.6451  -  .4237  +  4  -  4.5234  +  1.0216) 

=  .33995 


36  , .  9  ±  4  ,  „  5 


E[P?hl  = 


(R-A^)  i  B2-A1 
1  ‘  f A_-A . 1 f B*-B/V  1  n  B.-A. 


,  (19-13)  ,  22-13 

1  '  (17— 13)( 22-20)  In  20-13 


,  (19-17)  -  22-17 

+  (17-13)  (22-20)  1n  TtTTT 


1  -Th"!*!"’! 


=  1  -  (4.5)  (.2513144)  +  (.5)  (.5108256) 


=  .124498  . 


To  obtain  the  variance  of  Pkb,  the  following  formula  is  used 


Var[Pkb]  .  E[p£b]  -  (E[Pkb]}‘ 


=  .124498  -  ( . 33995) ^ 


=  .008932 


Let  a  =  /  Var[P  ]  *  .09451  . 
Kkb  KD 


Then  for  a  two-standard  deviation  confidence  interval  one  has 


pkb  =  Etpkb^  ±  2  °Pi 


=  .33995  +  (2)(. 09451) 


=  .33995  +  .18902. 


According  to  Chebyshev’s  inequality,  the  probability  is  at  least  75 
percent  that  Pth  will  lie  in  the  above  interval. 


%  .  *- 
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SECTION  III 


CONCLUSIONS 

A  methodology  has  been  developed  to  provide  confidence  Intervals  on  the 
probability  of  kill  due  to  blast,  P^,  in  the  absence  of  fragmentation  and 
aiming  error.  The  assumption  was  made  that  the  input  parameters  identifying 
were  subject  to  error  in  such  a  way  that  such  errors  could  be  described  as 
random  variables  defined  by  uniform  distributions.  Estimates  for  ECP^]  and 
VartPkfj]  were  provided  for  three  different  situations.  These,  in  turn,  were 
used  to  define  the  confidence  intervals  on  P^* 
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